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Abstract. We characterize the breaking of analyticity with respect to the replica 
number which occurs in random energy models via the complex zeros of the moment 
of the partition function. We perturbatively evaluate the zeros in the vicinity of the 
transition point based on an exact expression of the moment of the partition function 
utilizing the steepest descent method, and examine an asymptotic form of the locus of 
the zeros as the system size tends to infinity. The incident angle of this locus indicates 
that analyticity breaking is analogous to a phase transition of the second order. We 
also evaluate the number of zeros utilizing the argument principle of complex analysis. 
The actual number of zeros calculated numerically for systems of finite size agrees 
fairly well with the analytical results. 
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1. Introduction 

The replica method (RM) is a powerful tool in statistical mechanics for the analysis of 
disordered systems pLj . In general, the objective of RM is to evaluate the generating 
function 



for real replica numbers n G M (or the complex field C). Here, Z is the partition function, 
N represents the system size and (■ ■ ■) denotes the configurational average with respect 
to the external randomness which governs the objective system. Direct evaluation of eq. 
([T]) is generally difficult. However, eq. ([T]) can often be evaluated for natural numbers 
n = 1, 2, . . . G N in the thermodynamic limit as ^ oo. Therefore, in RM, one usually 
evaluates 



for n G N, in order to first obtain an analytic expression for (f){n) and then analytically 
continues this expression to n G M (or C). 

There are two known possible problems with this procedure. The first problem is 
multiple possible alternatives for analytic continuation [2J. Even if all values of (Z^) are 
provided for n E N, analytic continuation of (pNi^) from n G N to n G M (or C) is not 
uniquely determined, van Hemmen and Palmer conjectured that this may be the origin 
of the failure of the replica symmetric (RS) solution in the low temperature phase of 
the Sherrington-Kirkpatrick (SK) model [3]. The other issue is the possible breakdown 
of analyticity of (f){n). Although analyticity of (2'") with respect to n is generally 
guaranteed as long as N is finite, 0(n) = limTv^oo 0Ar('^) may fail to be analytic. This 
implies that if such a breaking of analyticity occurs at n = Uc < 1, then continuing the 
expression analytically from n G N to n G M will lead to an incorrect solution for n in 
the range of n < n^. 

In in [5], the authors developed an exact expression of the moment (Z^) for discrete 
versions of random energy models (DREMs) [6l [7]. The expression is valid for n G C 
and is useful for handling systems of finite size. Utilizing this expression, it can be 
shown that analyticity breaking of 0(n) actually occurs at a certain critical replica 
number n = Uc < 1 in the low temperature phase of DREMs. The uniqueness of the 
analytic continuation from n G N to n G C is guaranteed for DREMs. This means that 
the analyticity breaking of 0(n) with respect to n is the origin of the one step replica 
symmetry breaking (IRSB) which is observed for DREMs, and thus we are motivated 
to further explore its mathematical structure. 

This paper is written with this motivation in mind. Regarding log Z as the energy 
of the external randomness, the moment {Z"^) and generating function (j){n) are formally 
analogous to the partition function and free energy of non-random systems. This analogy 
leads us to characterize analyticity breaking in terms of the distribution of zeros of (Z") 
on the complex n plane, following the argument by Lee and Yang |8l [H [TOl [11] . Let us 
suppose that a partition function of a finite size system is expressed as a function of a 



0^(n) = llog (Z"), 



(1) 



0(n) = lim 0iv(n) 



(2) 
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certain control parameter (. As long as ( is real, the partition function does not vanish. 
However, there can exist zero points, which are generally termed the Lee- Yang zeros, 
throughout the complex ( plane. In general, the Lee- Yang zeros are apart from the real 
axis as long as the system size N is finite. However, when a certain phase transition 
occurs at a critical parameter value Cc G M, the distribution of the zeros becomes dense 
and approaches (c SiS N oo. The way in which this occurs characterizes the type of 
the phase transition. The main purpose of this paper is to apply such a description to 
the analyticity breaking of by evaluating the zeros of {Z") on the complex plane 
of the replica number n. 

This paper is organized as follows. In section Ej we briefly discuss the Lee- Yang 
zeros, utilizing a simple model. In section [3], which is the main part of this paper, we 
apply the Lee- Yang approach to examine the analyticity breaking of {Z"") of DREM. 
The exact expression for (Z") developed in [H [5] is utilized to perturbatively evaluate 
the zeros in the vicinity of the transition point employing the steepest descent method. 
This shows that the the distance of the closest zero to the real axis decays as 0{N~^^'^) 
when the system size is large. In the thermodynamic limit as oo, the incident 

angle of the locus to the real axis converges to 7r/4, indicating that analyticity breaking 
is analogous to a phase transition of the second order. However, it is also shown that 
the angle between the real axis and the line connecting the transition point with the 
kth zero from the real axis has a finite positive correction from it /A independently of A^ 
as long as k{= 1,2,.. .) is finite. The complex zeros are also numerically evaluated for 
several system sizes based on the expression, which is consistent with the asymptotic 
form to a reasonable precision. The final section is devoted to a summary. 

2. Complex zeros and analyticity 

2.1. Brief review of Lee- Yang zeros 

Partition functions of discrete systems of finite size are, in general, analytic with respect 
to their parameters because they are a summation of exponents of the Hamiltonian over 
all possible (but still a finite number of) states. In addition, they do not have zeros on 
the real axis for the same reason. 
A theorem by Weierstrass [T2] . 

"T/iere exists an entire function with arbitrarily prescribed zeros an provided 
that, in the case of infinitely many zeros, a„ ^ oo. Every entire function with 
these and no other zeros can be written in the form 



oo 




(3) 



where the product is taken over all an 7^ 0, the rUn are certain integers, and 
g{z) is an entire function." , 
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indicates that a partition function Z{() of a finite system can generally be expanded 
with respect to a parameter ( as 

oo ^ 

Z{C) = e^(^) - ^)e^+h&+-+^(^r- ^ (4) 

-I 0,n 
n=l 

where the are not on the real axis. Here, ( may be the inverse of the temperature /3, 
an external field, or any other parameters of the Hamiltonian. In most cases, partition 
functions are physically relevant only for real values of (. Eq. (jl]) yields an expression 
for the free energy. 



F(C) = --^logZ(C) 




where /3 > denotes the inverse temperature and C is a constant with respect to (, 
introduced so that the thermodynamic limit is well-defined. For systems of finite size, 
F{() is analytic in a neighborhood of the real axis since the a„ are located a finite 
distance away from the real axis. Therefore, systems of finite size do not exhibit phase 
transitions with respect to (. 

In the thermodynamic limit, however, the a„ can become dense and may approach 
the real axis. In such cases, the summation of the logarithm in eq. ([5]) is replaced by 
an integral with an integral contour which intersects the real axis. As a consequence, 
the free energy is expressed by a different analytic function on each side of the integral 
contour. This implies that a phase transition occurs at the intersection point as ( varies 
along the real axis. 



2.2. Simple example 

To illustrate the above scenario intuitively, we consider here a simple example. Let us 
suppose that the partition function of the example model can be written as, 

Z(C;a,Cc) = cosh^(C-Cc) (6) 
where a and (c are real positive numbers. The free energy is then 

i^(C;a,Cc) = -^logZ(C;a,Cc). (7) 

where 2a is introduced to make the thermodynamic limit of this model well-defined. 
Since coshz ^ gSign(Re(2))2^2 for \z\ 3> 1 holds , Z{(;aXc) behaves like 

ZiCa^OXc) 

^ ie^(C-C^)e(Re(C - Q) + le-^(^-^=)e(-Re(C - Q), (8) 

in the limit as a — > 0, where sign(a;) = x/\x\ for x 0, and B(x) = 1 for x > and 
vanishes for x < 0. Iie{z) denotes the real part of a complex number z. This shows that 
the free energy can be expressed as 

F(C; a ^ 0, Cc) = -^(C - Cc){e(Re(C - Cc)) - e(-Re(C - Cc))}, (9) 
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Figure 1. Zero points of the partition function Z{(;a,C,c) — cosh^(C — (c)- They 
become dense and approach the real axis in the limit as a — s- 



in this limit. This means that the analytic property of the free energy changes at the 
boundary Re(C) = Cc] more precisely, the first derivative of F{(] a, Cc) with respect to ( 
becomes discontinuous at Re(C) = Cc whereas the real part of the free energy remains 
continuous. If the parameter ( is the temperature or an external field, this indicates 
that a first order phase transition occurs at ( = Cc- 

This transition can be linked to the asymptotic behavior of zeros of ( in the complex 
plane, which arises in the "artificial thermodynamic limit" as a 0, as follows. Using 
the product expansion of coshz. 



k=l 

we can write 



cosh.=n(i+^^^4^) 



(10) 



^(C;«,Ce)^n(l + ^|^), (11) 

which implies that g{() and m„ are zero in the expression of eq.(jll). Eq. (|TT1) means that 
zeros of Z{(] a, Cc) do not occur on the real axis but are distributed on C = Cc=t (2/c — l)ai 
{k = 1,2, . . .) as is depicted in figure [H, where i = ^/—T. Therefore, the free energy of 
this model can be expressed as 

^(C;.,C)^-|g.og(i.^|^). (12, 

All the zeros C = Cc i {2k — l)ai correspond to branch points of the logarithm. 
Nevertheless, the free energy is still analytic along the real long as a is finite 

since there are no branch points on the real axis. 



Analytic properties of random energy models II 



6 



However, the situation changes in the hmit as a ^ since the distance of zeros 
from the real axis becomes infinitesimaL Actually, the summation has to be replaced 
by an integral thus: 



which implies that the integrand becomes singular at x = |C ~ Cc| when i{( — (c) or 
—i{C — Cc) is a non-negatie real number. As a consequence, this integral correctly 
reproduces eq. ([9]), which is singular along a line CI : C = Cc + "^i parameterized by 
V e M including the real critical parameter ( = (c- 

2.3. Remarks 

Two issues are noteworthy here. 

2.3.1. Argument principle and the number of zeros The first issue concerns the number 
of zeros inside a closed contour. For any complex function f[z) which is meromorphic 
inside a closed contour 7, the identity 



holds, where aj and hk are the zeros and poles inside 7, respectively. n{'y,p) represents 
the winding number of 7 around p. This formula is sometimes termed the argument 
principle [12]. Applying this to eqs. (jlj) and ([5]), one can evaluate the number of zeros in 
the region surrounded by 7 In most cases, analytically finding all zeros is difficult and 
one has to resort to numerical schemes such as Newton's method. A major difficulty of 
such approaches is to determine whether there has been a sufficient number of search 
trials. On the other hand, in some cases, the free energy density can be obtained 
in a computationally feasible manner in the thermodynamic limit. Therefore, one can 
estimate the asymptotic number of zeros by applying eq. f|T^ to the free energy density. 
This approach can be utilized to check whether or not sufficiently many zeros have been 
obtained. 

2.3.2. Incident angle and type of phase transition The second issue is the relation 
between the locus of the zeros and the type of the phase transition. In the simple model 
mentioned above, the locus that the zeros form in the limit as a ^ is the straight line 
CI : C = (c + vi (w e M), which is perpendicular to the real axis. This can be reproduced 
by dealing only with the limiting form of free energy 1^ as follows. Eq. (jnD means that 
in the limit as a ^ 0, the free energy is expressed by either of two analytic functions 
Fi{() = —{Tr/i3){( — (c) and -F2(C) = (^//5)(C ~ Cc), depending on the region to which C 
belongs. The critical condition for the selection is provided by Re (Fi(C)) = Re (F2(C)), 

I Eq. p4)) is analogous to Gauss's law for a two dimensional electromagnetic field. In this analogy, the 
free energy and zeros (poles) correspond to the electrostatic potential and point charges, respectively. 




(13) 




(14) 
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which yields the locus CI. Such an argument can be generalized to some extent. Let 
us suppose that the free energy in the thermodynamic limit can be expressed by either 
of two analytic solutions Fi(C) and F2(C) depending on (, and the first discontinuity 
appears at the level of the first derivative. This means that the two analytic solutions 
are expanded as 

FiiO = ~ {A, + A,{C - Ce) + A^iC - Cc? + ...), (15) 

and 

i^2(C) = [B, + i?i(C - Cc) + i?2(C - Cc)' + ...), (16) 

where Aq = Bq and Ai ^ Bi. A simple scenario implies that when the scale factor a is 
small but finite, the partition function can be asymptotically expressed as 

Z{C; a, Cc) ^ Pi exp + p2 exp 



fAo , (Ai + i?i)(C-Cc) A ^ ^ (Ai-i?i)(C-Cc) , ,\ 
oc exp ^- + J X cosh + ej , (17) 

for IC — Cc| a^^"^, where pi and p2 are prefactors subexponential with respect to a~^ 
and \9\ <^ a^^. Since exp {Ao/a + {Ai + i?i)(C — Cc)/(2a)) never vanishes, zeros in the 
vicinity of Cc, satisfying |C — Cc| ^ o^/^, come out from only the cosh part, and are 
expressed as Cfc ^ Cc ± (2A; - l)7ra/(Ai - Bi)i - 2ae/{Ai - Bi) (A; = 1, 2, . . .). In the 
limit of a — > +0, distances of contiguous zeros become infinitesimal and the locus of the 
zeros is parameterized as C — Cc + (f € K) as a|6'| — > 0. 

This argument indicates that the condition Re(Fi(C)) = Re(F2(C)) forms a locus 
of zeros in the thermodynamic limit and, as long as the transition is of the first order, 
which means that the first derivative of the free energy becomes discontinuous at a 
critical parameter value Cc, the locus makes an incident angle (defined in the upper-half 
plane hereafter) 7r/2 to the real axis. The angle, however, depends on the type of the 
transition [lOl [TTl [13] . Let us suppose a case of the second order phase transition, which 
corresponds to a situation Aq = Bq, Ai = Bi and A2 7^ B2 in the above setting. For 
small but finite a, an argument similar to the above yields an expression 

7.. .^ (Aq , A(C-Cc) , {A2 + B2){C-Cc? 

Z{C; a, Cc) oc exp — + + 

\ a a 2a 

V 2a 

for IC ^ Cc| ^ a^^^, implying that zeros in the vicinity of Cc are asymptotically expressed 
as C — Cc + '^exp (±(2/ — l)7r/4i) (/ = 1, 2, v E R) and the incident angle to the real axis 
is not 7r/2 but n/A and 37r/4. More generally, when the first discontinuity of the free 
energy comes out at the mth derivative, possible incident angles are limited to forms of 
{21 — 1)11 /{2m) {I = 1,2, ... ,m) [H]. In this way, the profile of the locus provides us 
with a useful clue for classifying the types of phase transition. 

Notice that the above argument crucially relies on an assumption that the 
subexponential prefactors pi and p2 do not vanish in the vicinity of Cc- In some cases. 
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however, either of them can vanish in the left or right side of Cc, which excludes some 
of the multiple branches of the incident angles. Later, we will see that this does occur 
in DREM. 



3. Complex zeros with respect to the replica number for a discrete random 
energy model 

Now, we are ready to employ the Lee- Yang type approach for characterizing analyticity 
breaking with respect to the replica number which occurs in DREM. 

3.1. Model definition and an exact expression of the moment of the partition function 

A DREM is defined by sampling 2^ energy states €1,62, ■■■ , independently from an 
identical distribution 

where and M = aN are positive integers, and Ei = i — M/2 for i = 0, 1, . . . , M. 

For a given realization of energy levels {ei}, the partition function of an inverse 
temperature P is defined as Z = ^^=1 exp (— /Jej). In [5], the authors showed that for 
Vn G C, VA^ and VM the moment of the partition function can be expressed as 



nM „ / M 



2N 



(Z") = / rfp(-p)-"-i 5^P(E.)e-> (20) 

r(-^) Jh ) 

where u = e~^, T{n) = — 2i sin(n7r)r(n) and the integration contour H is defined 
as shown in figure [2l T{z) is the Gamma function. From this expression, one can 
analytically determine the limit (f){n) = limjv^oo A^~^ log (2'"), which indicates the 
following behavior. Let us define a critical inverse temperature as 

\ log (1-/^2-^(1 -a-i))- log (/^2-^(l -«-!)), a >1, ^ ' 

where hl^^ijj) is the inverse function of the binary entropy h2{x) = — xlog2(x) — (1 — 
x) log2(l — x) for < X < 1/2. For /? < Pc, 0('^) is provided by either of 

0RS1 (n) = n (log 2 + a log (cosh ( f ) ) ) , (22) 



^'RS2(^) = ( log2 + «log ( cosh ( ^ ) ) ) , (23) 



or 

/ f f 

rAoon (n^ = (lop- 2 -I- ry lop" I cosh I ■ ^ 

depending on n. More precisely, there exists 3nRs > 1 such that (f){n) = 0rs2(^) for 
^ > ''^RS whereas 0(n) = 0rsi(^) for n < urs- On the other hand, the behavior is 
different for /? > P^- For n > = Pc/P, 0(^) = 0rs2(^) holds. However, 0(n) for 
n < Uc is described by neither 0rsi(^) nor 0rs2(^) but another solution 

, / X naP pc , . 

0irsb(^) = ^-tanh— . (24) 
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Figure 2. Integration contour to calculate Eg. ([20|) . 



Comparison with the rephca analysis indicates 0irsb(^) corresponds to the one 
step rephca symmetric (IRSB) solution [SIE!- For DREM, the uniqueness of analytical 
continuation from n G N to n G C is guaranteed by Carlson's theorem [15]. This 
indicates that the analyticity breaking a.t n = Uc for /5 > /5c is the origin of the IRSB 
transition. 

In addition to the advantage of yielding analytical expressions (1221) . (123|) and 
in the thermodynamic limit as ^ oo, eq. fl20l) is useful for numerically evaluating 
(Z^) for finite N because the necessary cost for the computation grows only linearly 
with A^. This property holds for Vn G C, which is advantageous in searching for zeros 
of {Z"") with respect to n G C in order to characterize transitions that occur in the limit 
as — > oo. 

3.2. Analytical results 

Let us apply the argument of the previous section to DREM. For simplicity, we focus 
here on the case of the IRSB transition aX n = n^. assuming a > 1 and (3 > (3c. 

3.2.1. Locus of zeros Eqs. ([231) and (El) indicate that (^V^n^) 0rs2(^) = 
(a/?V4) (l - tanh2(/5c/2)) > (^V^n^) 0irsb(^) = at n = = (3^/(3, whereas 
{d / dn) (f)^_S2in) = (5/(9?t,) 0irsb(^) = (a/3/2) tanh(/3c/2) holds. This implies that the 
IRSB transition is classified to the second order. The naive argument provided in the 
previous section means that the locus of zeros can be locally parameterized as 



where / = 1, 2, G R. However, this result must be corrected; the branch of / = 2, which 
corresponds to the incident angle 37r/4, never appears due to the following reason. 
For ^ 1, asymptotic evaluation of eq. fl20l) yields an expression 



n ~ ric + f exp (±(2/ — l)7r/4i) , 



(25) 
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derivation of which is shown in [5]. Evaluating dominant contributions in the integral 
by the steepest descent method [I6] indicates that the asymptotic expression is further 
simplified as 

(Zn) ^ / PRS2e^'^^s^(") +PiRSBe^*^^s«("), < Re(n) < 1, 

PiRSse^'^it^sBW^ 0<Re(n)<nc ^ ' 

depending on the position of n, where p^,s2 = (tt/ (2A^a;))^/^(cosh(?T,/3/2))~-'^ and Pirsb = 
{Na(3{nc — n))~^. This is because the absolute value of the integrand is maximized in 
the integral interval < y < h2^{l — a~^) for Uc < Re(n) < 1 while the right terminal 
point y = /i^^(l — Q^^) offers a unique dominant contribution for < Re(n) < ric- 
Equation (1271) means that, for ^ 1, there are no zeros of < Re(r;,) < Uc- Solving 
(Z^) = with respect to n perturbatively in the vicinity oi n = n^. for rtf. < Re(?7.) < 1 
under a condition of iV ^ 1 yields an asymptotic expression of the zeros 



~ ric + — ( 1 ) cosh ^'^ 



(3 \ 16k J \2 J M Na 



X exp (^±(^^ + ±- \og{8kn')^ i^ + o {N^'^') , (28) 



(fc = 1,2, . . .), where a (A^~^/^) represents contributions which are relatively negligible 
compared to N~^^'^. 

This result indicates that the distance from the real axis to the closest zero, rii, 
decays as 0(A^~^/^) and, in the thermodynamic limit N — ^ oo, the incident angle of the 
phase boundary converges to 7r/4. However, the angle between the real axis and the line 
connecting the transition point Uc with the kth zero, {k = 1,2, . . .), is larger than 
7r/4 by l/(8A;7r) log(8/c7r^) independently of as long as k is finite. This implies that 
for systems of finite sizes zeros in the vicinity of Uc are expected to be placed in the left 
side of the phase boundary. 

3.2.2. Number of zeros Eqs. fl23|) and flM|) can also be utilized to evaluate the 
asymptotic number of complex zeros. For this, let us consider a closed cycle 7 : Ci — >^ 
C2 — C3 — s> C4 as shown in figure [31 Application of the argument principle of the 
previous section to eqs. (j23l) and (1241) yields the asymptotic number of zeros inside 7. 7 
in figure [3] guarantees that the estimate does not vanish since the critical replica number 
= Pc/ P is in the range < n < 1. 
To simplify the analysis, we decompose the generating function as 

0;v(n) = ^ log (Z") = A(n) + iQ{n) (29) 

where A{n) and Q{n) are real functions. The argument principle indicates that 
N{2Tr)~^ J^dn{d/dn)fl{n), which is proportional to the total variation of Q{n) as n 
goes around 7, accords with the number of zeros inside 7. For A^ ^ 1, we substitute 
4>{n) into (j)i\f{n), which makes it possible to analytically evaluate the variation. 
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Figure 3. A cycle 7, the number of complex zeros inside which is evaluated utilizing 
the argument principle. 



(30) 



(31) 

(n) becomes 



Along Ci, both 0rs2(^) and 0irsb('^) lead to Q{n) = 0. Therefore, 

Ac,n = -n{0) = 0. 

Along C2, (piji) = (pRS2{n) = log(2) + a log (cosh {n(3 /2)) holds. This leads to 

Ac,fi = Im (0RS2 (1 + (27r//5)i)) - Im (0rs2 (1)) 
= Im (a log (— cosh(/9/2))) = avr, 

where Im(z) denotes the imaginary part of a complex number z. Along C3, 
real and therefore 

Acgfi = 0. 

Along C4, on the other hand, (pin) = 0irsb('^) holds, which yields 

Ac,n = Im (01RSB (0)) - Im (0irbs ((27r//5)i)) 

'a/3 ff3A2m\ f (3, 

T^^^H^J — J = -«7rtanh ^- 

These results indicate that the total number of zeros inside the cycle in the low 
temperature phase (3 > (3c can be asymptotically estimated as 



Im 



(32) 



(33) 



^est 



N aN 

— (Acifi + Ac.n + Ac,vi + Ac,n) 

ZTT 



(34) 



1 + e^- 
for > 1. 

This indicates that the number of the zeros inside 7 grows proportionally to A^. 
However, the total number of zeros over the complex plane of n is infinite because {Z^) 
is a transcendent function of n and has period 2t[ / (3 in the imaginary direction. This is 
in contrast to the case of the Lee- Yang zeros with respect to (3 for a typical sample of 
DREM [HllTj. In this case, the partition function Z{(3) becomes a polynomial of (3 and 
the number of zeros over the entire plane grows only linearly with respect to A^, which 
yields a continuous distribution of zeros after being averaged over {ej}. 
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Table 1. The number of complex zeros with respect to n in the cycle shown in figure 
|3]for a = 4 and /3 = 3/3c- The estimated number is calculated from eq. (|34p and the 
actual number v is evaluated by numerically searching for the zeros using eq. (j35p . 



3.3. Numerical results 

In order to justify the analytically obtained results for ^ 1, we numerically examined 
the zeros of (^"') = for several utilizing the expression of eq. ( l20l) . Unfortunately, as 
eq. (!20!) is transcendent with respect to n, solving for the zeros algebraically is difficult. 
Therefore, we resorted to an iterative numerical method to search for the zeros. 

We employed the secant method [17]. For solving an equation f{z) = 0, this scheme 
iterates the recurrence relation 

= fiz!)-Tiiy^''^^ (35) 
until \zt — Zt-i \ becomes smaller than a feasible discrepancy level, which was here chosen 
to be 10~^. This update rule is somewhat similar to that of Newton's method. Actually, 
the secant method can be regarded as an approximation of Newton's method, which is 
obtained by replacing {zt - Zt-i)/{f{zt) - f{zt-i)) with f'{zt) in eq. ([35]). A major 
advantage of this method is that there is no need to evaluate the first derivative. This 
property is useful when f{z) is complicated, which is the case for the current objective 
system 

To find all zeros inside 7 : Ci C2 ^ C3 ^ C4, we adopted the following strategy. 
We first covered the region inside 7 with a mesh of a fixed size in order to determine a 
set of initial values ni and n2, which are required when utilizing the secant method. Two 
adjacent points in the vertical direction on the mesh were taken as the initial points. 
In the iteration of eq. (!35ll . a calculated new point may step out of the region. In such 
cases, we changed the pair and tried the search again. Once we found a root ai, we 
replaced {Z"-) with {Z"') /{n — ai) and continued the calculation. After trying every pair 
of initial values, we reduced the mesh size and searched for the roots again. We finished 
these procedures when we were unable to find a new root. 

Table [H shows dependence of the number of zeros on the system size A^. This 
indicates that the number of zeros obtained numerically is reasonably consistent with 
the number obtained from analytical evaluation flM]) . implying that the number of the 
search trials has been sufficiently many. 

Figure H] shows results of the search for A^ = 5, 10, 15, 20 and 25. The curve 
represents the locus of the zeros in the limit a.s N —>■ 00, which was evaluated using 
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Figure 4. Complex zeros of the moment of the partition function of DREM for a = 4 
and P = SPc- The transition point in this case is determined to be Uc = (3c/ P =1/3. 
The curve represents the locus of zeros in the thermodynamic limit, which is evaluated 
using Re{(f>YiS2{n)) = R.e(0iRSB('^))- The incident angle of the curve to the real axis is 
7r/4, which indicates that the transition of 4>{n) is analogous to a phase transition of 
the second order. Numerical data (markers) approach the locus as the system size N 
increases from = 5 to 25 in steps of 5. 



the condition Re(0RS2('^)) = R'e(0iRSB(^))- As we expected, the incident angle of the 
locus to the real axis was 7r/4, indicating that the first discontinuity of <f){n) appears at 
the level of the second derivative. 

The markers represent data obtained from the numerical search. These show that 
the zeros approach the locus from the left side as N becomes larger, as we expected. 
Makers of figure [5] stand for zeros obtained from the asymptotic expression of eq. (127|) . 
These exhibit behavior qualitatively similar to that of the exact expression; but deviation 
is not negligible, which is presumably due to finite size effects. Nevertheless, the distance 
of Til to the real axis (Im(ni)) shows a fairly good consistency with the theoretical 
prediction of eq. (!28|) . validating our asymptotic analysis based on the steepest descent 
method (figure [6]). 

In conclusion, the overall consistency between the analytical and numerical 
evaluations justifies the current analysis based on eq. fl20|) . 

4. Summary 

In summary, we have characterized analyticity breaking with respect to the replica 
number n which arises in the discrete random energy model (DREM) by examining 
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1.8 




Re(n) 



Figure 5. Complex zeros obtained from the asymptotic expression of eq. for 
a = 4, /3 = 3/3c and TV = 5, 10, ... , 25. 



Exact 

Asymptotic 
Second Order 
First Order 




1.4 1.6 1.8 



2.2 2.4 2.6 2.8 



3.2 3.4 



log(N) 



Figure 6. Distance of ni to the real axis (Im(rii)) versus system size N. Markers, 
+ and X, represent the results for the exact expression of eq. ([20| and for the 
asymptotic one of eq. ([27| . respectively. The full line shows theoretical prediction 
Im(ni) = 5V7F/(4/3\/]VQ)cosh(/3c/2)sin(7r/4 + (1/8) log(87r2)), which is obtained by 
neglecting the o(N^^^^) contribution in eq. (|28p . while the correction (1/8) log(87r^) is 
omitted in the broken line. These support our theoretical predictions which indicate 
that the distance vanishes as 0{N~^^^) and positive corrections of the incident angle 
from 7r/4 exist for Uk independently of N as long as = 1, 2, ... is finite. 
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zeros of the moment of the partition function over the complex n plane. To do this, 
we utilized an exact expression of the moment of the partition function, which was 
introduced by the authors in a related paper [5]. The expression is valid for Vn G C 
and useful for numerically evaluating the moment of finite size systems. Taking the 
thermodynamic limit of the expression at sufficiently low temperatures indicates that 
analyticity breaking occurs at a critical replica number < < 1, which can be 
regarded as the origin of the one step RSB (IRSB) solution. Perturbatively evaluating 
the zeros in the vicinity of Uc based on the expression utilizing the steepest descent 
method shows that the distance from the real axis to the closest zero decays as 0{N~^^'^) 
when the system size is large. Examining the asymptotic form of the locus of the 
zeros in the thermodynamic limit N ^ oo implies that the transition is analogous to 
a phase transition of the second order. We also evaluated the asymptotic number of 
zeros inside a unit cycle shown in figure [3] based on the expression. Zeros numerically 
obtained for finite size systems are reasonably consistent with the analytical predictions. 

The approach developed here is also applicable to the standard continuous random 
energy model [18], and this approach yields results qualitatively the same as those for 
DREM. This implies that the IRSB transitions observed in a family of REMs can be 
generally characterized by the complex zeros in a similar manner. 

A natural question to ask is whether another type of RSB, full RSB (FRSB), can 
be characterized similarly by the complex zeros. Recently, one of the authors examined 
the zeros of tree systems in a vanishing temperature limit [19]. Although some of the 
systems are conjectured to exhibit FRSB at a certain critical number, numerical data 
about the zeros seem irrelevant to the FRSB transition. It is, therefore, desirable to 
explore other systems in order to assess whether or not the irrelevance of the zeros to 
FRSB is specific to tree systems. 
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